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ABSTRACT 

We investigate the production of electrons and positrons in the Milky Way within 
the context of dark matter annihilation. Upper limits on the relevant cross-section 
are obtained by combining observational data at different wavelengths (from Haslam, 
WMAP, and Fermi all-sky intensity maps) with recent measurements of the electron 
and positron spectra in the solar neighbourhood by PAMELA, Fermi, and HESS. We 
consider synchrotron emission in the radio and microwave bands, as well as inverse 
Compton scattering and final-state radiation at gamma-ray energies. For most values 
of the model parameters, the tightest constraints are imposed by the local positron 
spectrum and the final-state radiation from the central regions of the Galaxy. Accord- 
ing to our results, the dark matter annihilation cross-section into electron-positron 
pairs should not be higher than the canonical value for a thermal relic if the mass of 
the dark matter candidate is smaller than a few GeV. In addition, we also derive a 
stringent upper limit on the inner logarithmic slope a of the density profile of the Milky 
Way dark matter halo (a < 1.3 if m dm < 100 GeV and a < 1.8 if m dm < 10 TeV) 
assuming that {av) e ± = 3 x 10" 
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1 INTRODUCTION 



Dark matter can be indirectly detected through the 
signatures of standard mode l particles produced by its 
annih ilation or decay (see e.g. iBertone et al.ll2u05l ; iBertonel 
2010). A great deal of work has focused o n the emission of 
gamm a rays from the Ga lac tic ce ntr e (e.g. Berezinskv et al.1 



200£ 


Abdo et alj20ld; 


Acciari et alj|2010l: Essie et alJ200Sl. 


2010 


Ackermann et al. 


1201 ll: Abramowski et al. 20111). and 


galaxy clusters (e.g.lColafrancesco et alj2006t Jeltema et al.l 


2009; Ackermann et al 


2010a: Sanchez-Conde et all 2011: 


Pinzke et al.l l201ll). Prospects for indirect dark matter 



have also been 
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Morselli et al.1 2002 


; Ullio et al. 
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; Peirani et al. 2004 


; Prada et al. 
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; < 


Jesarini et al.1 12004 


Bergstrom et al. 2005a b; 


Profumo 


20051; Aharonian et al 


120061; Zahariias & Hooper 


20061; 


Bovarskv et al. 


2008; Pospelov et al.1 20081; 



detection in the microwave background 

considered by several authors (e.g. iBlasi et al.l 



200£ 


; iFornasa et al. 


20091; 


Bernal & Palomares-Ruiz 


201C 


; Abazaiian et al. 


2010 


; Cirelli et all 


201C ; 


Papucci & Strumia 201C 


1; Hooper & Goodenoughl 


2011; 


Hooper & Linden! l2011alll 


J; lAckermann et al.l |2012|. among 



Colafrancesc 
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2012) 



iPadmanabhan fc Finkbeine 



Lavalle 2010 



McQuinn 



Hiitsi et al 
Zaldarriagal 



2011 



wel l as X -ray (e.g 



20091 : iGalli et all 120091 ; ISlatver et al.1 120091; 



2011 



2003; 



2005: 



Mapelli et al.1 1 20061; IZhang et all 120061. 120071; ICholis et al 



Kanzaki et al 



Galli et al 



Dcla have et al 



Abazaiian et al. 2001 



200' 



Colafrancesco fc Meld l200ll; 

Borriello et al.1 



Bovarskv et~alj 120071. l200Sj ; IZavala et al.1 l201ll) radio 



many others), th e Milky Way satell it es (e.g. [ Baltz et all 
20001; iTvlerl |2002|; iBaltz fc Wa.il [iool; iHooper et al.1 1 2004 | 



Bcrgstro m fc Hooper 120061; ISanchez-Conde et al.l 120071; 
Strigari et al.ll2007l . 120081 ; IWood et al.ll2008l ; iMartinez et ail 



Aloisio et al 
[2009! 



2004; 



Bergstrom et al.l |2009|; iBqrriello et al.1 12009] ; llshiwata et al.l 
Forn engo et al ll 2l) multi-wavelength signatures 



(e.g. [Regis! | 2008|; iRegis fc Ulliol liooi ; IBertone et all |2009| 



Pato et al.l 120091 : iProfumo fc UIlicTbOlOl ; 
20101 ). 



Crocker et al.l 
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The recent results from indirect detection experiments 
in the solar neighbourhood have also suggested the possi- 
bility that such a signature has been seen. In particular, 
the PAMELA experiment has pointed a significant excess 
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of electrons and posit rons above the expect ed smooth astro- 
physical background (|Adriani et al.ll2009al ). If these results 
are interpreted in terms of dark matter annihilation, then 
an abundant population of high-energy e is being created 
everywhere in the Galactic dark matter halo, with the as- 
sociated final-state radiation (FSR), as well as synchrotron 
emission in the Galactic magnetic field and inverse Compton 
scattering (ICS) of the photons of the interstellar radiation 
field (ISRF). 

Although the currently most favoured explanation for 
the origin of Galactic positrons, traced b y the positron 
annih ilation emission line at 511 keV (see iPrantzos et al.l 
120111, for a recent review) is low-mass X-ray binaries 
I Weidenspointner et alj|2008h . and the local positron excess 
at high energies is mos t likely due to the contribution of 
nearby pulsars (see e.g. |Profumoll2012l ). several works have 
considered the possibility that dark matter annihilation 
makes a sizeable c ontribution to the p o sitron budget of the 
Milky Way (e.g. iBoehm et all 12004 iBoehm fc Ascasibarl 



Milky 
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Beacom et al 



Ascasibar et all 



2006 



2007 
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Picciotto fc Pospelovl |200 



IFinkbcincr 
Barger et al 



2006; iBeacom fc Yiiksell 120061; ISizun et all 
20071; 



Chen fc Takahashi 200 



Wcincr 
20091; 



9; 



Pospeli 



Bcrgstrom ct al 



Cholis et al 



2004 



2009|; lDonato et al.l2009l;lGrasso et aTll2009l; 
20od ; Mertsch fc Sarkarll2009l;lRegis fc Ullidl2009l;lYm et al 




Cirelli et al 



Malvshev et al 



20091; IChen et alj|20ld : iMeade et aLll201fj| ; ICline et al.ll201ll ; 



Vincent et al.ll2012l T 



The present work focuses on the astrophysical sig- 
natures of dark matter annihilation into electron-positron 
pairs, neglecting other processes, such as dark matter decay, 
or other annihilation products, such as protons and antipro- 
tons (whose contribution i s severely constraine d by recent 
observational data; see e.g. [Ad nam et al.ll2009bh . We try to 
impose robust, yet stringent constrains on the relevant cross- 
section by comparing the predictions of an analytic model of 
particle propagation with a multi-wavelength set of observa- 
tional data obtained from the literature. More precisely, we 
compare the expected emission from synchrotron radiation, 
ICS and FSR within the Milky Way with 18 maps of the sky 
at different frequencies: the Haslam radio map at 408 MHz, 
the 7-year data from the Wilkinson Microwave Anisotropy 
Probe (WMAP) in its 5 bands (23 GHz, 33 GHz, 41 GHz, 
61 GHz, and 94 GHz), and gamma-ray maps from the Fermi 
Large Area Telescope (LAT) binned in 12 different channels 
(from 0.3 to 300 GeV). A straightforward statistical criterion 
is used in order to mask the most obvious astrophysical sig- 
nals (i.e. emission from the galactic disc and point sources), 
and observational upper limits are derived from the remain- 
ing spherically-symmetric component. 

In addition to the photon data, we also consider the re- 
cent measurements of the local electron and positro n spec- 
tra p erformed by PAMELA jAdriani et all l2009al . l2010l . 
l201ll). Fermi jAckermann et all l2010bl . 120121 ). and HESS 
( Aharonian et al. 20081 ). As will be shown below, consider- 
ing the positron spectrum separately (rather than the com- 
bined electron+positron spectrum) yields a significant im- 
provement on the maximum value allowed for the positron 
injection rate or, equivalently, the dark matter annihilation 
cross- sect ion. 

Rather than focusing on a particular dark matter can- 
didate, we adopt a model-independent approach (see e.g. 



IWechakama fc Ascasibarl |201ll ). in which all the injected 
particles are created with the same initial energy Eo, of the 
order of the mass of the dark matter particle. Since this mass 
is usually much larger than the rest mass of the electron, 
electrons and positrons will be relativistic at the moment of 
their creation. However, they can efficiently lose their energy 
through different processes, such as ICS, synchrotron radi- 
ation, Coulomb collisions, bremsstrahlung, and ionization. 
Throughout this paper, we will often use the Lorentz factor 
7 to express the energy E — •ym e c 2 of the annihilation prod- 
ucts, where m c denotes the rest mass of electron, and c is 
the speed of light. We will first discuss the results obtained 
for a 'canonical' model of the Milky Way model and then 
explore the effects of varying the intensity of magnetic field, 
the diffusion coefficient, the ISRF, and the inner logarithmic 
slope of dark matter density profile. 

The remainder of this paper is structured as follows: 
Section [2] describes the procedure followed to estimate the 
electron-positron spectrum, the surface brightness profiles, 
and the parameters of the Milky Way model. Our analysis of 
the observational data is fully described in Section [3] (tables 
with precise numeric values are provided as an appendix), 
and Section [4] is devoted to the constrains on the dark mat- 
ter annihilation cross-section. The effect of the different as- 
trophysical parameters is discussed in Section [5] while Sec- 
tion[6]focuses on the constraints that one can impose on the 
slope of the dark matter density profile by assuming that 
dark matter particles are produced as thermal relics in the 
primordial universe. Finally, our main conclusions are suc- 
cinctly summarized in Section [7] 



2 MODEL PREDICTIONS 

2.1 Electron-positron propagation 

As in our previous work (|Wechakama fc Ascasibarl 1201 lh , 
the propagation of electrons and positrons through the in- 
terstellar medium (ISM) is determined by the diffusion-loss 
equation 



d An 
dtd^ 



&0> 7)^0 B > 7) 



d_ 

d-y 

Q(x yl ). 

We assume a diffusion coefficient of the form 
if (7) = Kol S 



+ 



(1) 



(2) 



independent of Galactic location. The values of Kq 
and S correspondin g to the three models discussed by 
Donato et all (|2004T ) are provided in Table Q] below. The en- 
ergy loss rate 



b(x, 7) 



d7 
~dt 



0,7) 



6i(x,7) 



(3) 



is a sum over the relevant physical processes, and the source 
term Q(x,y) represents the instantaneous electron-positron 
injection rate. 

Given e nough time (of the order of 100 Myr; c.f. 
Figure 2 in IWechakama fc Ascasibarl l201ll ) , the electron- 
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positron population will approach a steady-state distribu- 
(x, 7) = 0. Assuming that b(x, 7) varies smoothly 



tion 



' at dj 



in space, the particle spectrum fulfills the relation 

Q(x,l), (4) 



where 

y{x,i) = 6(7)^(35,7)- 



(5) 



Imposing ^{x, 7) = at infinity, one obtains the 
Green's function 



G(a;,7, a; s ,7 s ) = 



exp 



( \x-x s \ 2 \ 



(2ttAA 2 ) 



3/2 



e( 7 - 7s) 



(6) 



and either the image charges method or an expansion over 
the eigenfunctions of the linear differential operator may be 
used to derive t he Green's function for other boundary con - 
ditions (see e.g. iBaltz fc Edsiolll999l ; iDelahave et al.ll2009t ). 
The electron-positron spectrum is thus given by 



dn . . _ 1 
where the quantity 



d7 s / d^; 



exp 



|x-x s 

2AX A 



(2ttAA 2 ) 3/2 



AA 



A 2 (7) - A 2 ( 7 s) 



<?(xs,7s), (7) 



(8) 



is related to the characteristic diffusion length of the elec- 
trons and positrons, 7 S denotes their initial energy, and the 
variable A is defined as 



A 2 (7) 



2^(7) 
6(7) 



d7- 



(9) 



Considering the dark matter halo as a spherically- 
symmetric source, the spatial integral can be reduced to 
one dimension, and the electron-positron spectrum is finally 
given by the expression 

dn 1 ^(-2^) 

d7 (r '' 7) K7) (27rr 2 AA 2 ) 1/2 



d7 s / dr B r B exp 



exp 



VAA 2 



exp 



/ rr s \ 
V AAV 



2AA 2 , 

Q(r s , 7s )}. (10) 



2.2 Loss rates 



Electrons and positrons can lose their energy by several 
physical processes as they move through the ISM. We con- 
sider ICS of cosmic microwave background (CMB), starlight 
and infrared photons, synchrotron radiation, Coulomb col- 
lisions, bremsstrahlung, and ionization of neutral hydrogen 
atoms. 

The energy loss rates depend on the energy of the parti- 
cle. High-en ergy electrons and positrons mainly lose energy 
by ICS (e.g. ISarazinll 19991 ). The relevant loss function is 



6ics(7) = ^-^7 2 ^ad, 
3 m B c 



(11) 



where ctt is the Thomson cross section. The combined ra- 
diation energy density of the CMB, starlight (SL), and 



infrared (IR) light from thermal dust emission (see e.g. 
IPorter fc Strong) 120051 ; IPorter et all 120081 ) is represented by 
three grey bodies, 



4(JSB 



(12) 



where Ti and represent the effective temperature and 
the normalization of each component, respectively, and usb 
is the Stefan-Boltzmann constant. The cosmic microwave 
background is modeled as a perfect black body with tem- 
perature Tomb = 2.7 26 K (jFixsenl [20091 ) . and we follow 
ICirelli fc Pancil (|2009l ) for the two other components (see 
Tabled below). 

Synchrotron radiation is another important loss mech- 
anism at high energies. The expression for the loss rate is 
similar to that of ICS, substituting the radiation energy den- 
sity in equation (|ll|l by the energy density of the magnetic 
field, Ub = B 2 /(Stv), where B is the intensity of the mag- 
netic field: 



. . 4 OT 2 rr 

1(7) = o 7 U B - 

6 m e c 



(13) 



For lower-energy electrons and positrons, Coulomb in- 
teractions with the thermal plasma m ust be taken in to ac- 
count. The loss rate is approximately (R.ephacli 1979) 



6coui(7) » 1-2 x 10 2 n c 



1 + 



ln(7/n e 



75 



(14) 



where n c is the number density of thermal electrons. 

Collisions with thermal ions and electrons also 
produce radiation through bremsstrahlung. The loss 
rate due to bremsstrahlu ng can be approximated as 
|Blumenthal fc Gould|[l97rj| ) 



6brem(7) « 1-51 x 10~ 16 n c7 [ln( 7 ) + 0.36] 



(15) 



Additional energy losses come from the ionization of 
hydrogen atoms. The loss rate is given in lLongairl l|l98ll ). 



&ion(7) 



qt 



In 



7(7 2 - 1) 



r 



ln2+4 + i 

7 2 8 



1 - 



7 



(16) 



where wh is the number density of hydrogen atoms, q c is 
the electron charge, eo is the permittivity of free space, and 
I is the ionization energy of the hydrogen atom. The num- 
ber density of thermal electrons and neutral atoms can be 
expressed in terms of the total ISM gas density p g and the 
ionization fraction X; on clS 



and 



7lH = 



- ( 1 — ^ion ) , 



(17) 



(18) 



respectively. 



2.3 Source term 



Since the electrons and positrons in our model originate from 
the annihilation of dark matter particles, the instantaneous 
production rate at any given point can be expressed as 
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Q(r,j) = r) n dm (r) n dm «(r) (av) c 



dN„ 



-(7), 



(19) 



where ridm and ridm* denote the number densities of dark 
matter particles and anti-particles, respectively, {av) e ± is 
the thermal average of the annihilation cross-section times 
the dark matter relative velocity, and is the injection 

spectrum of electrons and positrons in the final state. For 
self-conjugate dark matter particles, ndm = ndm* — Pdm 
and rj — 1/2 in order to avoid double counting; else, 7id m = 
™dm- = |^ and 7? = 1. 

We assume that each annihilation event injects one 
electron and one positron with roughly the same energy 
70 ~ m dm /m c , 

dN, 



(1- 



-(7) = 2 5(7-70), 



(20) 



where 5(7 — 70) denotes a Dirac delta function. Although 
this is a rather coarse approximation, it has the advantage 
of being model-independent. For self-conjugate dark matter 
particles, we obtain 



Q(r,l) = 



Pdm(f) 



(av) e ± 5(7-70). 



(21) 



We consider a spherically-symmetric halo, described by 
a density profile of the form 



Pdm(f) = 



(22) 



where r s and p s denote a characteristic density and radius of 
the halo, respectively, and a is the inner logarithmic slope of 
the density profile. Local inhomogeneities that would boost 
the expected signal, such as small-scale dumpiness or the 
presence of subhaloes, are not taken into account. The shape 
of the dark matter density profile in the inner regions is 
far from being a settled question. N-body simulations sug- 
gest that, at least in the absence of baryons, the profile 
should be quite steep near the centre (a ~ 1), in apparent 
contradiction with observations. Traditionally, it has been 
argued that the presence of gas and stars makes the pro- 
file even steeper due to the effects of adiabatic contraction 
jBlumenthal et alj|l986l ). although some recent claims have 
also been made in the opposite direction (e.g. El-Zant et alj 
|2001| ; iMashchenko et al.M200rj ; lOh et al.1 l2010f ). Given "the 
current uncertainties, we have left the inner slope of the 
density profile as a free parameter of the model. 



2.4 Surface brightness profile 

Once the electron-positron spectrum is computed, the emis- 
sion coefficien10 for photons of frequency v is given by the 
integral 



e "( r > u } = ~^ f ^( r >7) i(7,^) d7, 



(23) 



of the electron-positron spectrum jp* (r, 7) times the specific 
luminosity v) emitted at frequency v by a single electron 
or positron with Lorentz factor 7. The intensity from any 

1 Energy radiated per unit volume per unit frequency per unit 
time per unit solid angle. 



given direction in the sky is simply the integral along the 
line of sight of the emission coefficient. Since we assume 
a spherically-symmetric source, it will only depend on the 
angular separation 8 with respect to the Galactic centre, 



e(r, v) ds, 



(24) 



where s represents the distance along the line of sight, and 
the radial distance r to the centre of the Milky Way at any 
point along the ray is 



r — \J x 2 + y 2 , 



(25) 
5 kpc (the 



with x = ss'mB, y = scosd — Rq, and Rq 
distance of the Sun from the Galactic centre). 

The contribution of synchrotron radiation, which dom- 
inates at low photon energies, can be estimated as (see e.g 
ISarazirj|i"999h 



i S p(7,^^%(7)], 



(26) 



where q c is the electron charge, B is the intensity of the 
magnetic field, and the function R(x) is defined as (e.g. 
iGhisellini et alJll98Sf ) 



R(x) = 2x 2 



K 



(X)*J (X)- §X -*|(X)} 



(27) 



In this expression, K refers to the modified Bessel function, 
and the normalized frequency 

X=A- (28) 



37 2 ^ c 

is expressed in terms of the cyclotron frequency 
q e B 

Vc = ~ • 

27rm e c 

At high photon energies (i.e. gamma rays), we consider 
the contributions of inverse Compton scat tering and final- 
state radiation. For ICS (e.g. ISarazi 



(29) 



jics(7i v) — 12-kot I J 
'0 



47 2 ?7 



with <tt denoting the Thomson cross-section and 
= 1 + 77 + 2r7ln?7 - 2r) 2 . 



(30) 



(31) 



Here, J is the mean intensity of the interstellar radiation 
field being scattered, which we represented as the sum of 
three grey bodies 

JW = ^ 



+ 



exp(/ii/ /fcTcMB) 



exp(hvo /kTs 



1 



exp(hvo/kTn 



1 



(32) 



where h is the Planck constant, and k is the Boltzmann 
constant. 

For FSR, the emission coefficient for photons of fre- 
quency v is given by 



e v {r,v) = 



hv 
4tt 



Pdm(r) 



d(av)i 



dv 



(33) 



with each annihilation event yielding a photon spectrum 
given by 
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Model 



K [cm 2 s" 1 ] 



MIN 

MED (canonical) 
MAX 



7.06 x 10 26 0.46 
1.67 x 10 25 0.70 
7.69 x 10 23 0.85 



Table 1. Three differen t mode ls of the diffusion coefficient pro- 
posed by iDonato et al ] j2004h . following the parameterization 

K(rt) = k 01 s . 



Model 



= 3481 K 



■An 
T m = 40.i 



K 



ISRF (I) 1.7 x 10- 11 
Canonical 2.7 X 10" 12 
ISRF (II) 8.9 x 10" 13 



7.0 x 10" 5 
7.0 x 10~ 5 
1.3 x 10~ 5 



a 


p s c 2 [GeV cm 3 ] 


r s [kpc] 


0.00 


2.346 


8.64 


0.20 


1.737 


9.56 


0.50 


1.042 


11.41 


0.70 


0.702 


13.08 


1.00 


0.349 


16.67 


1.20 


0.197 


20.33 


1.25 


0.169 


21.49 


1.50 


0.066 


29.81 


1.70 


0.025 


42.57 


1.90 


0.006 


70.30 



Table 3. Characteristic density and radius of the dark matter 
density profile fl'il'li as a function of its asymptotic logarithmic 
inner slope a. 



Table 2. Normalization of the grey-bod y models describing t he 
interstellar radiation field, adopted from lCirelli fc Pancil {2009). 



d(a"u}FSR a k 2 -2k + 2 

= (av) e ± In 

av iv v 



(1-k) 



(34) 



where a is the fine-structure cons tant and k = hu/rridr, 
(see e.g. IPeskin fc Schroededll995h . 



2.5 Astrophysical parameters 

The emission coefficient associated to final-state radiation 
is fully specified by the initial energy and injection rate of 
the electron-positron pairs, related to the nature of the dark 
matter particle (mass and cross-section) and the parameters 
describing the density profile of the Galactic halo. In con- 
trast, the photon intensity from the synchrotron and ICS 
emission also depends on the astrophysical parameters that 
determine the propagation and energy losses of the relativis- 
tic particles. We will first define a canonical model based on 
observations of the Milky Way and then investigate the effect 
of each individual component by varying the values of the 
adopted parameters. In all cases, we calculate the electron- 
positron spectrum as described in expression (|10|) , and then 
estimate the photon intensity according to expression (|24[) . 

Our canonical m odel assumes a dark matter density 
profile with a — 1 (jNavarro et al] Il997f ). r s — 17 kpc 
and p s c 2 = 0.35 GeV cm , co nsistent with dynamical 
models of the M ilky Way (e.g. iDehnen fc Binnevl 1 19981 : 
iKlvpin et al]|2002t ). The virial mass of the Galaxy is thus 
10 12 Mq, and the local dark matter density is Pdm(r©) c 2 = 
0.3 GeV cm -3 . The ISM is mainly composed of neutral hy- 
drogen atoms (Xj on = 0) with number density p g /m D ~ 
1 cm (IDehnen fc Binney|[l998l ; iFerrierdlioOll ; iRobin et al.l 
120031 ), and it is permeated by a tangled magne tic field whose 
intensity is B ~ 6 nG throughout t he Galaxy (jFerrierd l200ll : 
lBeckll200ll ; lAscasibar fc Dfaj|2010h . 

Apart from the canonical model, we consider the effect 
that the magnetic field, the diffusion coefficient, and the 
interstellar radiation field have on the synchrotron and ICS 
emission. The intensity of the magnetic field B is varied from 
1 to 100 i-iG. For the diffusion coefficien t (see equation [ 2l), we 
consider the three models discussed bv lDonato et al] ( 2004] ) 
and summarized in Table [T] We will also use three differ - 
ent models of the ISRF (adopted from lCirelli fc Panci|[2009l ) 
where the photon intensity is represented by three grey-body 



components (see equations 1121 and 132 . The normalizations 
and effective temperatures of the light emitted by the Galac- 
tic stars and dust are quoted in Table [2] 

Most importantly, we also investigate the effect of the 
inner slope of dark matter density profile on the production 
rate of electron-positron pairs. We vary the inner logarithmic 
slope a from 0.5 to 1.5. When varying 01 we also modify 
the characteristic density and radius in expression (|22|l so 
that the dark matter density at the solar radius is equal to 
0.3 GeV cm and the virial mass of the Galaxy is 10 12 M . 
The appropriate values of p s and r s are quoted in Table [3] 
for several values of the inner logarithmic slope a. 



3 OBSERVATIONAL DATA 

In order to constrain the production of relativistic electrons 
and positrons in the Milky Way, we consider observations of 
the whole sky at very different wavelengths. More precisely, 
the Haslam map in the radio band, the 5 WMAP channels 
at microwave wavelengths, and 12 energy bins of the Fermi 
LAT observations in the gamma-ray regime. The Haslam 
and WMAP maps are dominated by synchrotron emission, 
whereas Fermi traces ICS and FSR. 

The Haslam 408 MH z radio continuum all-sky map 
ijHaslam et al.lll98ll . fl982i ) combines data from four differ- 
ent surveys. The data were obtained from the archives of 
the NCSA ADIL in equatorial 1950 coordinates, and they 
were subsequently processed further in the Fourier domain 
to mitigate baseline striping and strong point sources. For 
the WMAP data, we take the full-resolution coadded tem- 
perature maps for each of the 5 frequency bands (23, 33, 
41, 6 1, and 94 GHz) cor responding to the 7-year observa- 
tions (jjarosik et al]|2011 ). T he Fer mi gamma-ray maps were 
computed bv iDobler etaD (|2010h from all "Class 3" (dif- 
fuse) photon events in the first-year data release. We use 
the 12 logarithmically-spaced frequency bands, from 0.3 to 
300 GeV, of the smoothed maps without point source sub- 
traction. 

Since we are interested in a spherically-symmetric com- 
ponent, we may follow a simple, conservative procedure in 
order to mask the emission from the Galactic disk and indi- 
vidual point sources without relying on any particular fore- 
ground model. For each frequency, we compute the average 
intensity 1(0) in 180 bins as a function of the angular sep- 
aration 6 from the Galactic centre. We also estimate the 
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Figure 1. Theoretical photon spectra of synchrotron radi- 
ation (dashed black lines), ICS (dotted blue lines) and FSR 
(solid magenta lines) for our canonical model with (av) e ± = 
3 X 10 -26 cm 3 s _1 and different injection energies, evaluated at 
10° from the Galactic centre. Grey bands illustrate the frequency 
ranges of Haslam, WMAP and Fermi. The observational data at 
9 = 10° are plotted as red triangles. 

standard deviation cr(8) within each bin, as well as the av- 
erage standard deviation 



where n = 180 is the total number of the bins. We then start 
an iterative procedure, where all pixels more than 3<r a ve 
away from 1(6) are discarded until convergence is achieved. 

This method seems to correctly identify and remove the 
most obvious structures in all but the two highest-energy 
Fermi bands, where the photon statistics is so poor that 
it is extremely difficult to distinguish diffuse emission from 
individual point sources. For these two bands, we opted to 
use the original average intensity Io(6) without applying any 
mask. Raw intensity maps, masked residual maps, i.e. / — 
1(6), and the average intensity 1(9) for each wavelength are 
shown in Figures IA11 IA21 and IA3I of Appendix [X] Numeric 
values of I (6), 1(6), and a (6) are quoted in Tables |AT1 lA"2l 
[All and El 

Besides these observational data, we also consider the 
energy spectra of cosmic-ray electrons and positrons in 
the solar neighbourhood; in particular, we use the com- 
bined electron+po s itron s pectrum measur ed by the Fermi 
jAckermann et al.1 l2010bl l and HESS (|Aharonian et al.1 
l2008r ) collaborations, as well a s the positron-only spec trum 
determined from Fermi data |Ackermann et aTTl2012l '). For 
P AMELA, we combine the electron-only spectrum obtained 
b vlAdriani et al.l J201ll) with the positron fraction discussed 
111 lid riani et alTl|2010] ) in order to derive the positron spec- 
trum. 



4 CONSTRAINTS ON THE DARK MATTER 
CROSS-SECTION 

Once the emission from the galactic disc and the most promi- 
nent point sources is excluded, the remaining spherically- 
averaged component can be used to place upper limits on 



Figure 2. Angular separation 9 that provides the upper lim- 
its for synchrotron, ICS and FSR emission. Lines and shadowed 
regions show the average and standard deviation across different 
wavelength channels, respectively. 



the cross-section for dark matter annihilation into electron- 
positron pairs. 

First of all, model intensities are computed according to 
the scheme described in Section [5] We consider the injection 
energy (i.e. the mass of the dark matter particle) as a free 
parameter and investigate values of the initial Lorentz factor 
70 between 2 x 10 3 and 2 x 10 7 , corresponding to injection 
energies Eo = 7om e c 2 from 1 GeV to 10 TeV. As an example, 
Figure [1] displays the results of our canonical Milky Way 
model for the synchrotron, inverse Compton, and final-state 
radiation contributions to the photon intensity at 10° from 
the Galactic centre, assuming a dark matter annihilation 
cross-section of {av} e ± — 3 x 10~ 26 cm 3 s . 

One can readily see that the Haslam radio map will be 
most sensitive to synchrotron emission by particles with an 
initial energy between 1 and 10 GeV, whereas WMAP data 
will cover the range Eo ~ 10— 100 GeV. On the other hand, 
the gamma rays observed by the Fermi LAT will constrain 
the maximum ICS and FSR emission allowed. The final-state 
radiation is sharply peaked at the injection energy, and it 
traces values of Eo between 1 GeV and 1 TeV. The inverse 
Compton spectrum is broader, and it features three distinct 
emission peaks, due to the scattering of CMB, starlight, and 
infrared photons. It is best suited to probe injection energies 
between 10 and 100 GeV, although the maximum intensity 
is always much smaller than the FSR peak, and therefore 
the associated constraints are going to be weaker. 

Since the value of the annihilation cross-section only 
sets the normalization of the spectra, and it does not alter 
its shape, it is relatively easy to set an upper limit by im- 
posing that the model intensities do not exceed the observed 
values (red triangles in Figure [TJ at any angular separa- 
tion 6. Not surprisingly, the tightest constraint will always 
be provided by a small value of 6, i.e. close to the Galac- 
tic centre. The dark matter density, and thus the injection 
rate, are higher there than anywhere else in the Galaxy. 
However, the observed intensity also reaches a maximum at 
6 = 0, and particles may diffuse from their injection point, 
effectively smoothing the density cusp. The angle that sets 
the maximum normalization that would be compatible with 
the observations (i.e. the upper limit of (av) e ±) is plotted 
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Figure 3. Surface brightness profiles of synchrotron, ICS, and FSR as a function of the angular separation 8 from the Galactic centre. 
Red triangles correspond to the mean observational intensity after discarding the contribution of the Galactic disk and prominent 
point sources as discussed in Section [3] Theoretical profiles are normalized to the maximum value of the annihilation cross-section (see 
Figure[4} allowed by these data. The angular separation that provides the tightest constraint - i.e. the tangent point between models and 
observations - is depicted in Figure [3] For synchrotron and ICS emission, the intensities obtained for the canonical Milky Way model 
are expressed in grey to black lines, where a darker colour represents a higher value of the injection energy Eg. The normalized intensity 
of FSR, shown as a green solid line, does not depend on Eq. 



in Figure [5] and the predicted surface brightness profiles of 
synchrotron, ICS, and FSR emission, normalized according 
to such prescription, are plotted in Figure [3] together with 
the observational data. 

Final-state radiation is produced at the very moment 
of pair creation, and thus it directly traces the positron in- 
jection profile, which is, in turn, proportional to the square 
of the dark matter density. Therefore, the intensity of the 
FSR emission does not depend on the injection energy of 
the particles or any astrophysical parameter other than the 
inner logarithmic slope a of the dark matter density profile. 
For this reason, the normalized surface brightness profiles of 
FSR depicted in Figure|3]do not depend on Eq. In our canon- 



ical model (where a = 1), and even more so if a > 1, the 
tightest constraints on the final-state radiation come from 
the very centre of the Galaxy (0 < 1°). 

For synchrotron and ICS emission, particle diffusion 
makes the intensity profile shallower, especially at high injec- 
tion energies. In general, one can say that photons of a given 
frequency trace electrons and positrons within a certain en- 
ergy range. If that range is close to Eq, these particles would 
have just been injected, and therefore the effects of particle 
propagation should be small, whereas, for large values of Eq, 
these electrons and positrons would have traveled a signif- 
icant distance from the point of injection, and the surface 
brightness profile will become considerably shallower. 
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This trend is indeed evident in Figure [3] surface bright- 
ness profiles become progressively shallower as one moves 
from Eo = 1 GeV to 10 TeV, and the effect is more pro- 
nounced for those channels that trace low-energy particles, 
i.e. Haslam, WMAP, and the lowest-energy Fermi bands. 
In the most extreme cases, diffusion keeps the electron- 
positron spectrum (and the ensuing intensity) roughly con- 
stant within the innermost 10 — 20° . For inverse Compton 
scattering, the tightest constraints on the annihilation cross- 
section come from 9 ~ 5° at all but the highest injection en- 
ergies, whereas for synchrotron emission the optimal angle 
increases steadily from 5 to 15° (see Figure |2J. 

In addition to the photons arriving from the cen- 
tre of the Milky Way, the dark matter annihilation cross- 
section {av) e ± is also strongly constrained by the observed 
abundance of relativistic electrons and positrons in the so- 
lar neighbourhood. In particular, we consider the recent 
measureme nts of the positron spec trum by the Fermi col- 
laboration (|Ackermann et al.1 [20121 ') and the PAMELA ex- 
periment. The latter data have not been published as 
such in the literature, but they can be trivially derived 
from the quoted positro n fraction jAdriani et all I2010T I 
and electron spectrum (|Adriani et aLT2011 ). Since the 
positron fraction is of the order of 10 percent or less at 
the energies below ~ 10 GeV, the constraints from the 
positron-only spectrum will be much tighter than those 
derived from the combined electron+positron data. For 
the sake of compari son, we also show these for PAM ELA 
dAdriani et al.ll201ll). Fermi dXckermann et al.ll2010bT ), and 



HESS ( Aharonian et alj|2008t ). Note that, in the latter case, 
the measurements are able to probe higher (~ TeV) energies, 
but it is not possible to discriminate between the electron 
and positron signatures. 

Our constraints are derived by imposing that the pre- 
dicted amount of electrons and/or positrons does not exceed 
the observed values for any Lorentz factor 7. Given the en- 
ergy dependence of the observed spectrum, [jg] . ~ E~ 3 , 
and the energy losses, b(E) ~ E 2 , the most restrictive con- 
straint comes from the spectrum near the injection energy, 
where propagation can be safely neglected and [gf|] , l ~ 
E~ 2 . The maximum production rate allowed by the 
data can then be expressed as 

dn ' 



Qo(r G ) < 6(70) 



dE 



(7o) 

obs 



and one arrives to the condition 



<<™)e±(7o) < 



Pdm(r ) 



6(70) 



dn 



dE 



(7o) 

obs 



(36) 



(37) 



in order not to overproduce the observed signal. 

The results are plotted in Figure 21 together with the 
upper limits on the dark matter annihilation cross-section 
derived from the comparison of the predicted synchrotron, 
ICS and FSR emission, assuming our canonical Milky Way 
model for particle propagation, with multi-wavelength ob- 
servations by Haslam, WMAP, and Fermi. As can be readily 
seen in the figure, the tightest constraints are provided by 
the final-state radiation from the Galactic centre for injec- 
tion energies above 20 — 30 GeV, whereas the positron spec- 
trum in the solar neighbourhood limits the production cross- 
section at lower energies. The constraints obtained from syn- 
chrotron and inverse Compton scattering are much weaker, 
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Figure 4. Upper limits on the dark matter annihilation cross- 
section derived by comparing the predicted synchrotron (red 
dashed line), ICS (blue dot-dashed line) and FSR (black solid 
line) emission with multi-wavelength observational data. The ar- 
eas shaded in blue show the constraints obtained from the mea- 
surements of the combined electron+positron spectrum at the 
solar neighbourhood by PAMELA, Fermi and HESS. The up- 
per limits obtained from the positron spectrum are shown by 
the red areas. The horizontal dotted line indicates the value 
(av) e ± = 3 X 1(T 26 cm s _1 . 



by about one order of magnitude, and similar to the up- 
per limits associated to the combined electron and positron 
spectrum in the solar neighbourhood. 

The typical value for thermal relics, {av) e ± = 3 x 
10 -26 cm s^ 1 , is ruled out for particle masses lighter than a 
few GeV. Both particle physics processes and astrophysical 
boost factors have previously been advocated to increase the 
current annihilation rate in the Milky Way by more than a 
factor of 10 with respect to the early universe. According to 
Figure 31 such models are excluded for any dark matter can- 
didate below the ~ TeV regime annihilating primarily into 
electron-positron pairs. Since our analysis involves a very 
conservative treatment of the astrophysical signal, merely 
excluding the emission from the disk and prominent point 
sources, it is expected that a deeper understanding of the 
astrophysical sources of electrons and positrons would make 
possible to probe the interesting region of the parameter 
space below (crv} e ± = 3 x 10~ 26 cm s -1 . 



5 EFFECT OF THE ASTROPHYSICAL 
PARAMETERS 

All the constraints represented in Figure [4] are based on the 
'canonical' Milky Way model discussed in Section [2.51 One 
should note that the most stringent results are provided by 
the final-state radiation from the Galactic centre and the lo- 
cal positron spectrum. These observables directly trace the 
instantaneous injection rate, and therefore they do not de- 
pend on the propagation parameters. However, the surface 
brightness profiles of synchrotron and ICS emission are sen- 
sitive to the precise values adopted for the intensity of the 
magnetic field, the diffusion coefficient, and the interstellar 
radiation field. The inner logarithmic slope of dark matter 
density profile has a very strong impact on the injection rate 
close to the centre, and thus it affects all the tracers con- 
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Figure 5. Upper limits on the dark matter annihilation cross- 
section from synchrotron and ICS, for different values of the mag- 
netic field (top), interstellar radiation field (middle), and diffusion 
coefficient (bottom) . Constraints from FSR and the local positron 
spectrum are indicated by the shaded area. 

sidered in the present work except the positron spectrum in 
the solar neighbourhood. 

Here we investigate the effect of the various astrophysi- 
cal parameters of our propagation model on the upper limits 
obtained for the dark matter annihilation cross-section. As 
we did for the canonical model, we consider different ini- 
tial energies Eo from 1 GeV to 10 TeV and compare the 
predicted emission with the full observational data set, but 
now we vary each of the astrophysical parameters in turn in 
order to assess their influence on the results. 

Let us start with the intensity of the magnetic field 



B. This parameter plays an important role in the energy 
losses, and it sets the total amount of energy that is radiated 
away as synchrotron emission. The top panel on Figure 
shows the upper limits derived by combining the constraints 
obtained from synchrotron and inverse Compton scattering. 
The results obtained B — 1, 6 (our canonical model), 30 and 
100 /iG are plotted as dotted, solid, dash-dotted, and dashed 
lines, respectively. All the other constraints (FSR and local 
positron spectrum) are independent of B, and are shown by 
the shaded area. 

Synchrotron constraints are most important at the low- 
est injection energies (Eo ~ 1 — 30 GeV), while the upper 
limits at higher initial energies (from ~ 30 GeV to 10 TeV) 
are due to ICS in the gamma-ray regime. The intensity of 
the magnetic field affects both processes in an opposite way: 
for low values of the magnetic field, all energy is lost by in- 
verse Compton scattering, and synchrotron emission is al- 
most irrelevant; as one increases the value of B, synchrotron 
constraints become more important at the expense of ICS 
emission. In the most extreme case (B = 100 /^G), gamma- 
ray constraints are negligible, and the upper limits derived 
from synchrotron radiation are well approximated by a pure 
power law. For any value of the magnetic field, the upper 
limits derived from FSR and the positron spectrum are al- 
ways more stringent than those derived from synchrotron 
and ICS. 

In the middle panel of Figure O we investigate the up- 
per limits of synchrotron and ICS for different models of the 
interstellar radiation field. As mentioned in Section [2. 5 1 wc 
adopt ed the parameterization proposed by ICirelli fc Pancil 
(2009j) in terms of three black-body components. The tem- 
peratures and normalizations of each component are sum- 
marized in Table [21 The effect of the ISRF is similar to that 
of the magnetic field, but in the opposite direction: a higher 
photon density results in a larger amount of energy being 
lost by inverse Compton scattering rather than synchrotron 
emission. Nevertheless, for reasonable values of the model 
parameters, the upper limits on (av) do not vary by more 
than a factor of three. 

As shown on the bottom panel of Figure [5] the ef- 
fect of the diffusion coefficient is even smaller. The upper 
limits are slightly more stringent when the electrons and 
positrons are allowed to travel a shorter distance from the 
place where they were injected, but the diffe rence between 
the th ree propagation models discussed by iDonato et al.l 
(2004]) is barely noticeable. 

Thus, we conclude that our results are not severely af- 
fected by the astrophysical uncertainties associated to par- 
ticle propagation. In particular, synchrotron and inverse 
Compton constraints are weaker than the upper limits from 
final-state radiation and the local positron spectrum for any 
value of the magnetic field, the interstellar radiation field, 
or the diffusion coefficient. 

This conclusion is also valid for any value of the inner 
slope q of the dark matter density profile, but this parameter 
plays a very important role in setting the actual constraints 
on (av). We have investigated several values in the interval 
< 01 < 2 (the appropriate values of p s and r s are quoted in 
Table[3} , and we report in Figure[6]the upper limits obtained 
from the comparison of the predicted final-state radiation, 
synchrotron, and inverse Compton scattering emission for 
a = 0.50, 1.00, 1.25, and 1.50 with our multi-wavelength 
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Figure 6. Upper limits on the dark matter annihilation cross- 
section for different values of the inner logarithmic slope a of the 
dark matter density profile. On the top panel, the constraints are 
derived from the angle that provides the most stringent limit (see 
Figurc[2]for the canonical case a = 1), whereas all the constraints 
on the bottom panel are obtained from the observed emission at 
8 = 10° from the Galactic centre. In both cases, black and red 
lines represent the limits associated to FSR and Synchrotron+ICS 
emission, respectively. The constraints from the local positron 
spectrum (independent on a) are shown by the shadowed areas. 



observational data set. Results for FSR and the combina- 
tion of synchrotron and ICS emission are plotted separately. 
Constraints from the local positron spectrum are indepen- 
dent of alpha and are shown as a solid area. 

The top panel of the figure shows the upper limits 
obtained by the same procedure applied to the canonical 
model, i.e. choosing the angular separation 9 that provides 
the tightest constraint. Not surprisingly, larger values of 
a result in lower values of 9. The constraints from FSR 
and Synchrotron+ICS emission come from innermost 1° for 
a >0.5 and a >1.25, respectively. 

Since the particle production rate near the centre of 
the Milky Way increases dramatically with the value of the 
inner slope of the density profile, this is, by far, the most 
relevant astrophysical parameter. For a > 1.25, the upper 
limits from final-state radiation are always stronger than 
the constraints provided by the local positron spectrum, and 
they rule out a cross-section larger than 3 x 1CF 26 cm 3 s _1 
for any dark matter candidate lighter than ~ 100 Gev. On 
the contrary, if the dark matter density profile of the Milky 
Way was shallow, with a logarithmic slope a < 0.50, the 
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Figure 7. Upper limits on the inner logarithmic slope of dark 
matter density profile a, obtained by imposing that FSR from 
the Galactic centre does not overproduce the observed signal (ac- 
cording to the Fermi data) for a thermal dark matter relic (i.e. 
(crv) e ± = 3 X 10" 26 cm 3 s" 1 ). 

positron spectrum in the solar neighbourhood would provide 
the most stringent limits on dark matter annihilation. 

One may remove the dependency of the results on the 
precise shape of the dark matter density profile by fixing 9 = 
10° when comparing model predictions with observational 
data. As shown in the bottom panel of Figure El we find, 
in agreement with previous w ork (e.g. ISerpico fc Zahariias! 
120081 : lAckermann et a.1.1 120121 ) . that the uncertainty associ- 
ated to the precise value of a reduces to about a factor of 2 
when the comparison is restricted to the photon intensity at 
9 — 10° . While this is therefore a good choice when the goal 
is to provide a conservative upper limit on the dark matter 
annihilation cross-section, we would like to stress that any 
prior knowledge of the dark matter density profile may lead 
to much stronger constraints if the inner slope was steeper 
than a = 1, as evidenced in the upper panel. 



6 CONSTRAINTS ON THE INNER SLOPE OF 
THE DENSITY PROFILE 

As pointed out in I Ascasibar et ail (|2006l ). the photons from 
the central region of the Galaxy contain information on both 
the dark matter annihilation cross-section and the shape of 
the density profile. By assuming a given value of the cross- 
section, one can constrain the value of a from the total inten- 
sity and the morphology of the observed surface brightness. 

In this work, we will focus only on the total intensity 
from the central bin (i.e. one degree around the Galactic cen- 
tre) in order to derive a robust upper limit. More detailed 
constraints could be obtained from the shape of the surface 
brightness profiles at different wavelengths once the astro- 
physical contribution is adequately subtracted. We set the 
dark matter annihilation cross-section into electron-positron 
pairs to the value expected for a thermally-produced relic, 
{av) e ± — 3 x 10 -26 cm 3 s _1 , and compute the value of a 
for which the predicted final-state radiation rises above the 
observed level. 

The corresponding upper limit is plotted in Figure [7] as 
a function of the initial energy Eo associated to the mass of 
dark matter candidate. Our results show that, for a thermal 
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relic with m^m < 100 GeV, the dark matter density profile of 
the Milky Way must be shallower than a ~ 1.3 in order not 
to overproduce the observed signal. It is worth noting that, 
since FSR only depends on the injection rate, our constraint 
on the inner logarithmic slope a is independent on the other 
astrophysical parameters. Moreover, the fact that we are 
considering the total gamma-ray emission, without taking 
into account the contribution of astrophysical origin, implies 
that this is a conservative upper limit, and therefore we can 
conclude that, if dark matter particles annihilate primarily 
into electrons and positrons (or, more generally, any lepton 
pair), any scenario where the Milky way features a steep 
density profile (due to e.g. adiabatic contraction) may be 
firmly ruled out. 



7 SUMMARY AND CONCLUSIONS 

We have investigated the constrains on the dark matter an- 
nihilation cross-section into electron-positron pairs by com- 
paring the predictions of an analytic model of particle propa- 
gation with a multi- wavelength set of observational data ob- 
tained from the literature. We have compared the expected 
emission from synchrotron radiation, inverse Compton scat- 
tering and final-state radiation within the Milky Way with 
18 maps of the sky at different frequencies: the Haslam radio 
map at 408 MHz, the 7-year data from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) in its 5 bands (23 GHz, 
33 GHz, 41 GHz, 61 GHz, and 94 GHz), and gamma-ray 
maps from the Fermi Large Area Telescope (LAT) binned 
in 12 different channels (from 0.3 to 300 GeV). A straight- 
forward statistical criterion has been followed in order to 
mask the most obvious astrophysical signals (i.e. the emis- 
sion from the galactic disc and prominent point sources), 
and observational upper limits are derived from the remain- 
ing spherically-symmetric component. In addition, we have 
also imposed that the predicted abundance of electrons and 
positrons in the solar neighbourhood does not exceed the 
measurements by PAMELA, HESS and Fermi. Our main 
results can be summarized as follows: 

(i) The constraints from synchrotron and inverse Comp- 
ton emission are always weaker than those from final-state 
radiation from the Galactic centre and the positron spec- 
trum in the solar neighbourhood. This result is valid for all 
the values of the magnetic field, diffusion coefficient, and 
models of the interstellar radiation field that we have con- 
sidered. 

(ii) If the density profile of the Milky Way halo is steep 
(a > 1), the annihilation cross-section is most tightly con- 
strained by final-state radiation, and the exact value of 
inner logarithmic slope plays a crucial role in the upper 
limit on {<rv) e ±. This scenario is hardly consistent with 
any dark matter candidate lighter than ~ 100 GeV and 
(av) e ± = 3 x 10~ 26 cm 3 s" 1 . 

(iii) If the density profile of the Milky Way halo is rela- 
tively shallow (a ^ 1), the upper limit on the cross-section 
is set by the local positron spectrum for low values of the 
injection energy. Combining both types of messenger (pho- 
tons and positrons) is thus of the utmost importance in this 
case. Considering the positron spectrum separately makes 
possible to rule out cross-sections above 3 x 10 -26 cm 3 s _1 
for dark matter particles lighter than a few GeV. 



Let us conclude by noting that the current upper limits 
are close to - or have just reached - the expected anni- 
hilation cross-section for a thermal relic. A better under- 
standing of the production of positrons and gamma rays by 
astrophysical sources should thus lead to the detection of an 
indirect signal from dark matter annihilation into electrons 
and positrons, providing at the same time an exquisite probe 
of the distribution of dark matter in the innermost regions 
of the Milky Way halo. Otherwise, the most straightforward 
version of the leptophillic dark matter scenario would be 
ruled out completely. 
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Figure Al. Haslam and WMAP intensity maps 1(1, b) in Galactic coordinates (left), masked residual maps 1(1, b) — 1(8) (middle), and 
spherically-averaged intensities (right). Dotted blue lines represent the original mean intensity Io(S), while solid red lines correspond to 
the final intensity 1(9) after discarding the outliers (black areas in the masked residual maps). 
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Figure A2. Fermi intensity maps 1(1, b) in Galactic coordinates (left), masked residual maps 1(1, b) — 1(9) (middle), and spherically- 
averaged intensities (right). Dotted blue lines represent the original mean intensity Iq(9), while solid red lines correspond to the final 
intensity 1(9) after discarding the outliers (black areas in the masked residual maps). 
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averaged intensities (right). Dotted blue lines represent the original mean intensity Iq(9), while solid red lines correspond to the final 
intensity 1(9) after discarding the outliers (black areas in the masked residual maps). For the last two bands, we opted to use the original 
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1.9 


6.7 


60.5 


9.7 


8.7 


2.2 


1.2 


2.0 


1.6 


3.1 


1.7 


2.3 


3.1 


1.7 


3.2 


3.3 


2.1 


7.1 


70.5 


8.7 


7.6 


1.6 


1.7 


1.7 


1.3 


1.0 


1,5 


2.1 


3.8 


1,5 


3.1 


3.9 


1.8 


6.8 


80.5 


9.0 


7.5 


1.7 


6.8 


1.5 


1.4 


5.8 


1.2 


2.0 


5.3 


1.1 


3.0 


4.9 


1.0 


6.6 


90.5 


8.0 


7.2 


1.9 


4.1 


1.4 


1.3 


3,1 


1.3 


2.2 


3.0 


1,1 


3.3 


3.0 


1.8 


7.0 


100.5 


7.6 


7.0 


1.7 


3.5 


1.8 


1.6 


3.0 


1.8 


2.3 


3.0 


1.9 


3.1 


3,5 


2.8 


6.9 


110.5 


7.5 


6.7 


1.7 


3.7 


1.5 


1.6 


3.1 


1.1 


2,5 


2.9 


1.6 


3.6 


3.2 


2.3 


7,5 


120.5 


7.2 


6.5 


1.8 


2.6 


1.4 


1.5 


2.0 


1.2 


2,1 


1.8 


1.2 


3,1 


1.8 


1.4 


7,1 


130.5 


7.5 


6.9 


2.1 


2.6 


1.8 


1.8 


2.0 


1,5 


2.7 


1.7 


1,1 


3.7 


1.6 


1.4 


7.7 


140.5 


8.3 


7.7 


2.3 


1.1 


2.8 


2.2 


3.1 


2.2 


2.7 


2.7 


2.0 


3.7 


2.7 


2.2 


8.1 


150.5 


9.4 


9.0 


2.0 


5.3 


3.0 


2.1 


1.1 


2.1 


2.8 


3.6 


2,5 


1.0 


3,5 


2.7 


8.2 


155.5 


9.9 


9.8 


2.2 


5.5 


3.8 


2.1 


1.1 


3.2 


2.9 


3.6 


2.8 


3.8 


3. 1 


3.0 


8.0 


160.5 


10.6 


10.5 


1.7 


6.1 


4.8 


2.5 


4.8 


4.0 


2.7 


4.2 


3.9 


3.6 


4.2 


4.1 


7.3 


162.5 


10.7 


10.7 


1.6 


6.1 


5.3 


2.5 


4.6 


4.3 


2.7 


4.0 


3.9 


3,1 


4.1 


4.1 


7.1 


164.5 


10.7 


10.7 


1.5 


6.2 


5.6 


2.6 


1.7 


1.6 


2.8 


1.1 


1.1 


3.1 


1,5 


1,5 


7.3 


166.5 


10.9 


10.9 


1.6 


6.8 


6.0 


2.8 


5.3 


1.8 


2.8 


1.7 


1. 1 


3,5 


5.3 


5.2 


7.1 


168.5 


11.2 


11.2 


1.7 


7.1 


5.9 


2.3 


5.2 


4.8 


2.9 


1,1 


4.2 


3.6 


4.6 


4.5 


7,1 


170.5 


12.1 


12.1 


1.8 


10.6 


6.2 


1.9 


8.5 


4.9 


2.9 


7.7 


4.9 


4.1 


8.1 


6.2 


8,5 


171.5 


12.4 


12.4 


1.8 


10.5 


7.0 


2.3 


8.0 


5.7 


3.0 


7.1 


5.1 


1.1 


7,5 


7.0 


8,5 


172.5 


12.7 


12.7 


2.0 


24.3 


7.8 


2.2 


21.5 


5.8 


2.3 


19.9 


5.1 


3,5 


18.6 


6.2 


7.9 


173.5 


12.6 


12.6 


1.9 


14.5 


9.0 


2.4 


10.8 


6.6 


2,5 


9.1 


5,5 


3.3 


8.2 


6.0 


7.6 


174.5 


12.3 


12.3 


1.9 


10.4 


9.9 


2.7 


7.6 


7,1 


2.6 


6,1 


6.3 


3.2 


6.7 


6.7 


6.7 


175.5 


12.0 


12.0 


1.6 


9.3 


9.2 


2.5 


6.6 


6.6 


2,5 


5.3 


5.3 


3.0 


1.9 


1.9 


6.1 


176.5 


12.0 


12.0 


1.4 


10.1 


9.8 


3.0 


7.1 


7.1 


2.8 


5.8 


5.8 


3.0 


5.1 


5.1 


5.9 


177.5 


12.4 


12.4 


1.3 


12.5 


14.7 


2.3 


9.7 


9.3 


3.8 


8.3 


8.2 


4.1 


9.0 


9.0 


7.7 


178.5 


12.4 


12.4 


1.2 


12.4 


14.7 


2.5 


8.9 


8.9 


3,5 


7,1 


7,1 


3,5 


6.9 


6.9 


6.2 


179.5 


12.7 


12.7 


0.9 


12.4 


12.4 


2.9 


8.9 


8.9 


2.7 


7.0 


7.0 


2.6 


6.0 


6.0 


1.9 



Table Al. Haslam and WMAP: Original observational mean intensity Io(8), final mean intensity after discarding the outliers 1(9), and 
standard deviation u(9), in units of GeV cm -2 s — 1 Hz -1 sr _1 . 
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WMAP 




Fermi 








94 GHz 


0.3-0.5 GeV 


0.5-0.9 GeV 


0.9-1.7 GeV 


1.7-3.0 GeV 


8 


I o (0) 1(8) a(9) 


I (9) 1(8) a(9) 


I (d) 1(9) a(9) 


I (9) 1(9) a(8) 


I (9) 1(9) a(9) 


[°] 


xl0~ 17 


xl0~ 30 


xl0~ 30 


xl0~ 30 


xl0~ 30 



22.2 1500.2 1506.8 34.3 989.9 1016.7 18.7 713.6 725.8 9.1 367.7 371.1 4.6 

24.5 1295.6 1353.7 37.3 774.3 790.1 16.6 447.3 462.2 9.8 230.1 245.2 4.7 

21.4 1037.0 992.3 37.0 553.9 508.8 19.8 275.9 219.7 9.9 139.7 116.5 4.9 
21.2 826.8 689.0 34.2 418.4 304.9 16.7 211.3 129.4 8.8 106.7 66.6 4.0 
21.8 680.6 460.1 36.1 340.2 208.4 17.8 174.1 107.2 9.1 86.8 53.3 4.3 

19.2 585.4 358.1 34.3 293.0 170.4 16.8 146.5 82.5 8.3 71.8 36.7 4.4 

18.3 509.1 289.2 32.1 258.2 133.7 16.1 129.9 56.2 9.3 63.4 26.9 4.1 

17.8 470.1 213.9 38.6 243.3 101.3 19.3 126.2 47.8 7.7 61.7 23.8 3.9 

18.1 415.0 184.3 35.3 211.3 90.0 17.5 107.7 42.6 8.0 52.1 21.5 3.8 

18.6 367.1 165.0 33.0 181.9 80.7 16.7 89.5 38.4 8.1 42.9 19.5 3.6 

19.0 341.4 150.7 31.9 167.5 74.2 17.0 83.5 36.8 8.6 39.8 18.3 3.9 

17.5 285.3 132.9 30.0 137.5 61.5 14.8 67.8 30.3 7.7 32.2 14.8 3.5 
18.5 271.2 123.1 31.7 130.1 56.9 15.8 63.7 28.2 8.0 31.5 13.5 3.5 

18.5 272.5 117.7 35.1 135.2 55.0 16.8 70.1 27.2 8.3 34.5 13.2 3.9 

19.4 264.6 110.3 34.7 132.4 51.7 16.7 66.6 25.7 8.7 32.2 12.3 4.0 

20.5 241.1 106.6 33.4 118.3 49.9 16.6 57.9 23.7 8.1 28.0 11.5 4.0 
18.4 193.1 89.0 27.0 93.5 40.1 12.4 46.3 19.7 6.7 22.7 9.6 3.2 
17.4 160.3 76.2 21.4 78.9 36.0 11.1 39.1 17.7 6.1 18.6 8.3 2.7 

15.9 115.6 60.7 19.6 56.0 27.9 10.0 27.8 13.2 5.1 13.4 6.4 2.4 

15.2 87.4 50.2 16.8 41.5 22.6 8.3 19.8 10.6 4.5 9.7 5.2 2.2 

16.1 65.4 45.1 18.4 30.1 20.0 9.3 14.4 9.5 4.8 6.8 4.5 2.3 
15.1 63.3 43.4 16.7 29.2 19.1 8.4 13.6 8.7 4.3 6.5 4.1 2.0 
14.9 63.3 39.2 16.0 29.2 17.1 7.7 14.0 8.0 3.9 6.5 3.7 1.8 

15.6 59.4 37.4 15.7 27.4 16.3 8.2 13.0 7.8 4.6 6.0 3.6 2.1 

15.4 57.8 37.4 15.6 25.9 16.6 7.9 12.0 7.9 4.3 5.7 3.8 2.1 

16.6 54.1 36.8 15.0 25.2 15.9 7.6 11.9 7.4 4.0 5.6 3.5 2.0 
16.6 51.3 36.3 13.8 23.7 16.0 7.0 11.4 7.4 3.5 5.5 3.6 1.9 

17.6 57.2 40.2 15.6 26.1 17.5 7.7 12.3 8.2 3.8 5.8 3.8 1.8 
18.9 70.5 49.4 20.4 32.8 21.4 9.6 15.7 10.2 4.9 7.1 4.9 2.4 

19.0 83.2 60.0 26.1 39.3 27.2 12.5 19.1 13.4 6.6 8.7 6.1 2.9 

18.8 89.8 71.3 29.5 42.7 31.9 14.1 20.3 14.7 6.9 9.2 6.9 3.3 

16.9 100.4 84.8 31.2 46.4 39.5 15.4 21.6 18.6 7.7 9.8 8.5 3.4 
16.8 126.4 94.1 36.9 63.9 41.8 16.9 29.0 20.8 8.9 13.7 9.2 3.9 

16.8 156.4 102.5 38.7 92.6 47.2 19.3 62.2 22.1 9.2 33.5 9.8 3.9 

17.3 142.8 109.6 36.7 68.4 49.1 18.0 29.7 23.2 9.1 13.9 10.8 4.1 
17.3 141.5 121.4 34.7 65.7 51.8 14.7 30.2 25.2 7.7 13.7 11.4 3.7 

19.9 173.4 129.9 33.0 83.6 60.1 16.0 41.2 28.7 7.9 19.3 13.5 3.7 

19.5 194.1 138.3 31.3 96.7 62.7 14.3 47.6 30.7 7.5 22.6 14.1 3.6 

18.1 209.0 154.7 33.0 107.2 68.0 13.2 57.2 35.8 8.2 27.6 15.2 3.4 
17.5 212.6 164.2 29.3 106.9 79.6 16.3 55.3 39.7 7.8 26.1 17.3 3.9 

15.7 203.9 172.7 26.6 95.7 82.0 13.2 44.9 41.2 7.6 20.5 19.0 3.7 
14.7 196.7 179.6 25.3 90.4 85.5 12.7 42.5 42.3 7.6 19.4 19.4 3.6 
13.9 196.4 193.0 29.5 92.9 92.9 14.6 45.1 45.5 8.7 20.7 20.5 4.1 
18.5 203.1 203.1 26.9 98.7 98.7 14.9 49.5 53.3 9.2 23.0 23.1 4.5 
14.9 211.8 211.8 21.3 103.2 103.2 13.2 52.3 55.7 9.3 24.9 24.9 4.0 

13.2 217.6 217.6 10.8 105.3 105.3 7.0 54.5 54.5 6.4 26.2 26.2 2.1 



0.5 


1097.8 


1132.3 


1.5 


287.4 


219.6 


2.5 


140.3 


66.9 


3.5 


114.6 


52.3 


1.5 


85.1 


39.9 


5.5 


69.9 


25.3 


6.5 


86.1 


20.6 


7.5 


66.6 


20.8 


8.5 


61.9 


18.0 


9.5 


43.2 


15.7 


10.5 


41.6 


12.5 


12.5 


33.5 


9.1 


14.5 


39.5 


11.7 


16.5 


36.6 


14.3 


18.5 


39.0 


15.0 


20.5 


27.8 


13.2 


25.5 


22.9 


11.1 


30.5 


20.9 


7.0 


40.5 


9.3 


3.5 


50.5 


9,1 


1,1 


60.5 


7.3 


4.9 


70.5 


7.3 


1.0 


80.5 


7.7 


2.2 


90.5 


6.1 


3.9 


100.5 


7.7 


6.2 


110.5 


6.9 


5.1 


120.5 


4.5 


3.5 


130.5 


4.1 


3.7 


140.5 


6.8 


5.6 


150.5 


8.6 


7.0 


155.5 


8.3 


7.9 


160.5 


10.3 


10.1 


162.5 


10.9 


10.8 


164.5 


12.3 


12.3 


166.5 


14.6 


14.4 


168.5 


13.5 


13.3 


170.5 


20.0 


17.8 


171.5 


19.7 


19.7 


172.5 


27.8 


17.7 


173.5 


20.9 


19.0 


174.5 


20.4 


20.5 


175.5 


16.4 


16.4 


176.5 


18.1 


18.1 


177.5 


26.8 


26.6 


178.5 


22.4 


22.4 


179.5 


21.0 


21.0 



Table A2. WMAP and Fermi: Original observational mean intensity Iq(0), final mean intensity after discarding the outliers 1(8), and 
standard deviation cr(8), in units of GeV cm~ 2 s _1 Hz -1 sr — 1 . 
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Fermi 

3.0-5.3 GeV 5.3-9.5 GeV 9.5-16.9 GeV 16.9-30.0 GeV 



6 

[°] 


He) 


1(8) 
xlO -31 


a(8) 




1(8) 
xlO -31 


a(8) 


Io(<>) 


1(8) 
xl0~ 31 


a(8) 


h(9) 


1(8) 
xlO -31 


a(8) 


0.5 


1554.0 


1574.6 


19.2 


482.4 


488.9 


7.5 


188.8 


195.6 


3.1 


74.4 


75.7 


1.8 


1.5 


950.3 


941.4 


17.8 


311.6 


321.6 


7.5 


118.7 


116.7 


3.1 


45.9 


44.7 


2.2 


2.5 


577.1 


468.4 


20.9 


201.6 


163.2 


7.3 


71.7 


60.3 


3.7 


29.3 


23.3 


1.9 


3.5 


443.4 


278.5 


17.8 


158.6 


101.2 


7.1 


51.5 


37.5 


3.4 


23.3 


19.1 


1.7 


4.5 


360.0 


220.5 


18.4 


129.0 


82.1 


7.9 


42.7 


27.8 


3.6 


19.3 


16.0 


2.0 


5.5 


305.0 


181.7 


17.6 


106.7 


60.1 


7.3 


38.6 


16.9 


3.2 


17.3 


11.9 


1.9 


6.5 


274.4 


149.0 


15.7 


93.0 


53.1 


6.1 


36.6 


15.6 


3.4 


15.6 


9.4 


2.1 


7.5 


263.3 


109.0 


21.7 


89.9 


47.3 


6.7 


36.2 


15.5 


3.3 


14.9 


9.1 


2.0 


8.5 


221.0 


94.7 


17.7 


78.7 


41.4 


8.0 


30.5 


13.9 


3. 1 


13.4 


7.2 


2.2 


9.5 


182.7 


85.8 


15.3 


67.2 


36.6 


8.2 


24.7 


11.9 


3.1 


11.3 


6.1 


1.9 


10.5 


169.4 


78.0 


17.2 


63.0 


33.2 


7.8 


22.4 


11.2 


3.2 


10.4 


5.7 


1.8 


12.5 


140.5 


66.4 


16.2 


52.6 


27.1 


6.7 


20.3 


10.0 


3.2 


8.8 


1.8 


1.8 


14.5 


132.9 


65.9 


16.5 


49.9 


25.6 


6.6 


19.5 


9.5 


3.1 


9.1 


5. 1 


1.9 


16.5 


146.3 


61.7 


16.4 


52.4 


24.0 


7.1 


18.9 


8.6 


2.5 


8.8 


1.6 


1.7 


18.5 


136.3 


53.9 


16.5 


50.5 


22.6 


6.6 


18.0 


8.5 


2.8 


7.5 


3.9 


1.7 


20.5 


116.8 


51.4 


17.0 


44.6 


20.4 


5.9 


17.0 


8.5 


3.0 


7.9 


1.0 


1.7 


25.5 


94.5 


42.6 


15.2 


35.2 


17.1 


5.3 


12.9 


6.0 


2.5 


6.0 


2.8 


1.5 


30.5 


79.0 


36.8 


10.6 


29.7 


14.5 


5.0 


10.9 


6.3 


2.2 


5.2 


3.0 


1.4 


40.5 


57.8 


29.6 


10.6 


21.2 


11.2 


4.7 


8.1 


4.8 


2.2 


3.6 


2.2 


1.1 


50.5 


39.8 


21.8 


8.7 


15.7 


9.5 


1.3 


5.9 


4.1 


2.0 


2.7 


1.9 


1.2 


60.5 


28.3 


19.9 


9.5 


11.6 


8.2 


3.6 


1.5 


3.3 


1.7 


2.1 


1.6 


1.1 


70.5 


27.0 


18.1 


9.0 


11.1 


7.8 


1.1 


1.6 


3.3 


1.9 


2.1 


1.5 


1.1 


80.5 


27.7 


16.5 


7.5 


11.3 


7.2 


3.8 


4.5 


3.0 


1.7 


2.2 


1.7 


1.2 


90.5 


25.7 


16.0 


8.6 


9.7 


6.9 


3.5 


3.8 


2.9 


1.8 


1.8 


1.5 


1.0 


100.5 


23.6 


16.2 


8.6 


9.2 


6.7 


1.0 


3.6 


3.0 


1.9 


1.7 


1. 1 


1.1 


110.5 


23.8 


15.7 


8.4 


9.3 


6.9 


3.7 


3.9 


3.0 


1.9 


2.1 


1.7 


1.2 


120.5 


22.6 


15.4 


7.9 


8.8 


6.5 


3.4 


3.5 


2.7 


1.7 


1.7 


1.5 


1.0 


130.5 


24.3 


17.8 


9.2 


9.9 


7.8 


4.3 


4.0 


3.0 


1.9 


1.9 


1.8 


1.2 


140.5 


30.0 


21.3 


11.2 


10.9 


8.6 


1.5 


1. 1 


3.7 


2.1 


2.0 


1.7 


1.2 


150.5 


35.2 


26.1 


12.5 


12.9 


10.7 


5.3 


1.6 


3.8 


2.1 


1.9 


1.7 


1.2 


155.5 


38.6 


29.0 


14.1 


13.9 


12.0 


5.8 


1.7 


1. 1 


2.2 


2.3 


2.3 


1.3 


160.5 


39.4 


35.3 


15.1 


13.6 


12.5 


4.7 


5.4 


5.3 


2.4 


2.7 


2.6 


1.5 


162.5 


54.9 


39.4 


18.1 


17.6 


14.2 


6.0 


6.2 


5.8 


2.4 


2.7 


2.4 


1.4 


164.5 


134.4 


42.9 


17.0 


39.3 


16.1 


5.8 


8.8 


5.6 


2.7 


2.9 


2.4 


1.5 


166.5 


55.7 


46.1 


18.1 


18.6 


15.5 


5.1 


6.9 


6.3 


2.3 


2.9 


2.8 


1.6 


168.5 


54.8 


46.9 


16.9 


19.9 


17.6 


5.6 


7.0 


6.6 


2.5 


3.3 


3.0 


1.6 


170.5 


82.7 


53.7 


14.5 


30.4 


21.6 


6.2 


11.2 


6.8 


2.7 


5.3 


3.2 


1.4 


171.5 


94.7 


54.1 


14.7 


35.3 


22.6 


6.2 


13.0 


6.7 


2.7 


5.9 


3.3 


1.5 


172.5 


111.4 


58.4 


17.4 


41.5 


22.9 


6.1 


16.2 


7.4 


2.9 


7.7 


3.6 


1.6 


173.5 


103.5 


65.2 


15.4 


37.0 


22.5 


5.9 


15.0 


8.7 


2.8 


7.2 


3.7 


1.8 


174.5 


79.0 


69.2 


11.6 


27.6 


24.0 


6.5 


10.2 


9.4 


2.1 


1.6 


3.9 


1.8 


175.5 


75.4 


74.7 


14.8 


26.1 


24.8 


5.6 


9.2 


9.0 


2.2 


3.9 


3.9 


1.3 


176.5 


81.8 


81.8 


17.2 


28.4 


28.3 


6.1 


9.9 


9.2 


2.1 


1.3 


1.3 


1.5 


177.5 


93.7 


98.6 


18.4 


32.7 


31.9 


7.2 


11.2 


11.0 


2.4 


4.8 


3.9 


1.4 


178.5 


105.6 


109.0 


22.0 


32.5 


32.4 


7.1 


10.3 


10.3 


1.5 


4.6 


4.2 


1.7 


179.5 


116.3 


116.3 


14.9 


29.1 


29.1 


3.6 


8.6 


8.6 


0.5 


3.8 


3.8 


0.8 



Table A3. Fermi: Original observational mean intensity Io(8), final mean intensity after discarding the outliers 1(8), and standard 
deviation cr(8), in units of GeV cm~ 2 s~ x Hz -1 sr -1 . 
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Fermi 

30.0-53.3 GeV 53.3-94.9 GeV 94.9-168.7 GeV 168.7-300.0 GeV 



e 




1(9) 


a(6) 


HO) 


1(6) 


a(B) 


h(8) 


Io(0) 


[°] 




xlO -32 






xlO -32 




xl0~ 32 


xlO -32 


0.5 


253.5 


263.1 


13.8 


130.5 


134.0 


8.7 


75.5 


56.9 


1.5 


189.2 


183.6 


14.5 


96.3 


101.4 


9. 1 


50.7 


31.4 


2.5 


144.4 


113.8 


12.9 


63.8 


50.9 


8.6 


32.1 


13.4 


3.5 


115.9 


76.8 


12.3 


AAA 


37.2 


9. 1 


22.3 


5.4 


1.5 


97.8 


58.0 


13.2 


37.4 


33.1 


8.4 


17.5 


4.8 


5.5 


86.2 


47.6 


11.4 


34.1 


27.0 


7.5 


12.4 


5.6 


6.5 


71.5 


39.8 


9.9 


32.1 


15.3 


9.3 


9.6 


5.7 


7.5 


71.0 


32.1 


12.8 


32.5 


17.3 


9.1 


10.0 


7.5 


8.5 


64.3 


29.7 


12.5 


30.1 


18.1 


9.3 


9.5 


6.6 


9.5 


55.4 


27.6 


11.1 


28.2 


14.3 


8.6 


9.6 


5.2 


10.5 


55.1 


28.4 


11.9 


29.8 


14.0 


9.0 


11.9 


1.7 


12.5 


48.5 


25.3 


11.8 


27.8 


12.6 


8.4 


10.7 


4.6 


14.5 


44.0 


20.9 


11.4 


20.1 


10.1 


6.9 


8.8 


1.9 


16.5 


43.7 


21.0 


10.2 


17.2 


7.3 


6.0 


8.2 


1.2 


18.5 


39.7 


21.9 


11.1 


17.9 


8.0 


7.0 


9.1 


3.7 


20.5 


39.2 


22.7 


11.4 


13.4 


7.8 


6.1 


9.5 


2.5 


25.5 


28.9 


16.5 


9.8 


14.4 


9.1 


6.6 


7.1 


3.5 


30.5 


21.8 


13.7 


9. 1 


10.9 


6.3 


6.0 


6.3 


2.0 


40.5 


19.1 


13.4 


8.6 


8.5 


6.6 


6.2 


1.2 


2.0 


50.5 


11.1 


8.0 


6.7 


6.8 


4. 1 


5.0 


4.6 


1.1 


60.5 


9.4 


7.7 


6.9 


6.2 


5.0 


5.5 


2.4 


0.8 


70.5 


11.4 


8.9 


8.2 


5.9 


1.8 


5.1 


3.6 


2.3 


80.5 


11.6 


8.6 


7.2 


5.9 


1.1 


5.1 


2.8 


1.1 


90.5 


8.9 


7.9 


7.7 


5.5 


4.8 


5.5 


2.4 


1.5 


100.5 


8.3 


7.5 


6.5 


5.2 


4.8 


5.4 


2.7 


0.9 


110.5 


9. 1 


8.3 


7.9 


5.1 


1.0 


5.2 


2.7 


1.1 


120.5 


8.3 


7.7 


7.0 


5.6 


5.3 


5. 1 


2.6 


1.2 


130.5 


11.5 


10.6 


8.5 


5.0 


4.5 


5.4 


2.4 


1.4 


140.5 


9.6 


8.7 


8.0 


6.0 


4.9 


5.4 


3.5 


1.5 


150.5 


10.9 


9.3 


7.9 


6.1 


5.2 


5.4 


2.8 


1.1 


155.5 


9.6 


8.4 


7.5 


5.3 


5.1 


5. 1 


3.3 


1.2 


160.5 


11.3 


10.6 


8.3 


5.7 


5.2 


5.7 


3.1 


2.7 


162.5 


9.8 


9.0 


7.6 


6.0 


5.0 


5.9 


3.4 


2.6 


164.5 


12.4 


11.6 


9.1 


6.9 


6.3 


6.0 


3.8 


1.5 


166.5 


12.8 


11.6 


7.5 


5.6 


5.5 


5.5 


3. 1 


1.7 


168.5 


11.6 


11.5 


8.7 


6.6 


5.9 


5.8 


3.5 


1.8 


170.5 


20.9 


13.3 


9.5 


12.1 


5.9 


5.2 


5.3 


2.6 


171.5 


23.4 


12.9 


9.1 


13.7 


6.3 


5.1 


6.9 


3.4 


172.5 


27.9 


12.0 


9.1 


18.2 


8.0 


6.3 


11.0 


1.7 


173.5 


26.4 


14.0 


8.5 


16.7 


7.9 


6.8 


9.6 


1.9 


174.5 


19.3 


17.7 


9.7 


9.1 


8.0 


6.7 


5.5 


3.3 


175.5 


17.6 


16.8 


9.6 


6.3 


6.3 


5.5 


6.5 


1.3 


176.5 


16.3 


16.3 


9.5 


5.6 


5.6 


5.3 


7.3 


0.9 


177.5 


17.8 


17.8 


9.1 


6.1 


6.4 


5.8 


1.6 


3.5 


178.5 


25.0 


25.0 


10.3 


6.9 


6.9 


6.0 


1.3 


5.2 


179.5 


29.7 


29.7 


9.3 


6.3 


6.3 


4.0 


0.2 


3.4 



Table A4. Fermi: Original observational mean intensity Io(9), final mean intensity after discarding the outliers 1(0), and standard 
deviation cr(9), in units of GeV cm~ 2 s~ x Hz -1 sr -1 . 
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